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Abstract 

Power corrections to the Q 2 behaviour of the low-order moments of both the longitudi- 
nal and transverse structure functions of proton and deuteron have been investigated 
using available phenomenological fits of existing data in the Q 2 range between 1 and 
20 (GeV/c) 2 . The Natchmann definition of the moments has been adopted for dis- 
entangling properly target-mass and dynamical higher-twist effects in the data. The 
leading twist has been treated at next-to-leading order in the strong coupling constant 
and the effects of higher orders of the perturbative series have been estimated using 
a renormalon-inspired model. The contributions of (target-dependent) multiparton 
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correlations to both 1/Q 2 and l/Q A power terms have been determined in the trans- 
verse channel, while the longitudinal one appears to be consistent with a pure infrared 
renormalon picture in the whole Q 2 -range between 1 and 20 (GeV/c) 2 . Finally, the 
extracted twist-2 contribution in the deuteron turns out to be compatible with the 
hypothesis of an enhanced ci-quark parton distribution at large x. 



1 Introduction 



The experimental investigation of deep-inelastic lepton-hadron scattering has provided a 
wealth of information on the occurrence of Bjorken scaling and its violations, giving a decisive 
support to the rise of the parton model and to the idea of asymptotic freedom. Quantum 
Chromodynamics (QCD) has been thereby proposed as the theory describing the logarithmic 
violations to scaling in the asymptotic region and its predictions at leading (LO) and next-to- 
leading (NLO) orders have been nicely confirmed by the experiments. However, in the pre- 
asymptotic region the full dependence of the hadron structure functions on the squared four- 
momentum transfer, Q 2 = q-q, is affected also by power corrections, which can originate from 
non-perturbative physics and are well beyond the predictive power of perturbative QCD. 
An important tool for the theoretical investigation of the Q 2 behaviour of the structure 
functions is the Operator Product Expansion (OPE): the logarithmic scale dependence is 
provided by the so-called leading twist operators, which in the parton language are one- 
body operators whose matrix elements yield the contribution of the individual partons to 
the structure functions. On the contrary, power corrections are associated to higher-twist 
operators which measure the relevance of correlations among partons (see, e.g., [|TJ). 

In case of unpolarised inelastic electron scattering the nucleon response is described 
by two independent quantities: the transverse F^(x,Q 2 ) and the longitudinal F^(x,Q 2 ) 
structure functions, where x = Q 2 /2Mu is the Bjorken variable, with M and v being the 
nucleon mass and the energy transfer in the target rest frame. Systematic measurement of 
the transverse structure function of the nucleon, F^(x,Q 2 ), (more precisely, of the proton 
and the deuteron [Q, ]3|, fi], ||, |(|) have been carried out in the kinematical range 10~ 4 < x < 1 
and for Q 2 values up to several hundreds of (GeV/c) 2 ; thus, various phenomenological fits of 
the data sets are presently available, like those of Refs. 0, |7|, ||. As for the longitudinal to 
transverse (L/T) cross section ratio, R^ T (x, Q 2 ) = <7l(x, Q 2 )/a^(x, Q 2 ), the experimentally 
investigated kinematical range is 0.0045 < x < 0.7 and 1 < Q 2 (GeV/c) 2 < 70 [|, |, g 
£|; however, the sparse and still fluctuating data are not yet sufficient to put significative 
constraints on the phenomenological fits (see Refs. [|H], |11|), particularly in the low-x region 
(x<0.3). 

Various analyses of power-suppressed terms in the world data on both FJ?(x, Q 2 ) 
and F^(x,Q 2 ) structure functions have been already carried out; they are based either on 
the choice of a phenomenological ansatz [fT2| , T3| or on renormalon-inspired models [H, IT5J, 
adopting for the leading twist the LO or NLO approximations. Very recently, also the 
effects of the next-to-next-to-leading order (NNLO) have been investigated on F^(x,Q 2 ) 
and R^ T (x,Q 2 ) [0 as well as on the parity-violating structure function xF^(x,Q 2 ) [pT7| , 



measured in the neutrino and antineutrino scattering off iron by the CCFR collaboration 
||18|| . Present analyses at NNLO seem to indicate that power corrections in F^(x,Q 2 ), 
R^/ T (x, Q 2 ) and xF^(x, Q 2 ) can be quite small. However, in such analyses the highest value 
of x is typically limited at ~ 0.75 and also the adopted Q 2 range is quite restricted (Q 2 > 
5 -=-6 (GeV/c) 2 at x ~ 0.75) in order to avoid the nucleon-resonance region; therefore, in such 
kinematical conditions the power corrections represent always a small fraction of the leading- 
twist term (corrected for target-mass effects). Furthermore, when power corrections are 
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investigated in a narrow Q 2 -window, it is clear that variations of the logarithmic dependence 
of the leading twist can easily simulate a small power-like term as well as a significative 
sensitivity to the higher-twist anomalous dimensions cannot be easily achieved. 

We point out that the smallness of the sum of the various dynamical higher-twists 
does not imply the smallness of its individual terms and therefore the question whether 
present data are compatible with an alternate-sign twist expansion is still open. To answer 
this question, we extend in this paper the range of values of Q 2 down to Q 2 ~ 1 (GeV/c) 2 
in order to enhance the sensitivity to power-like terms. For such low values of Q 2 we expect 
that the effects of target-dependent higher twists (i.e., multiparton correlations) should show 
up, thanks also to the contributions of the nucleon- resonance region and the nucleon elastic 
peak. The inclusion of the nucleon-resonance region is clearly worthwhile also because of 
parton-hadron duality arguments (see, e.g., Ref. [PHI ). 

In this paper we analyse the power corrections to the Q 2 behaviour of the low- 
order moments of both the longitudinal and transverse structure functions of proton and 
deuteron using available phenomenological fits of existing data in the Q 2 range between 
1 and 20 (GeV/c) 2 , including also the SLAC proton data of Ref. || covering the region 
beyond x ~ 0.7 up to x ~ 0.98. The Natchmann definition of the moments is adopted 
for disentangling properly target-mass and dynamical higher-twist effects in the data. The 
leading twist is treated at NLO in the strong coupling constant and, as far as the trans- 
verse channel is concerned, it is extracted simultaneously with the higher-twist terms from 
the data. The effects of higher orders of the perturbative series are estimated adopting the 
infrared {IB) renormalon model of Ref. ||14j| , containing both 1/Q 2 and 1/Q 4 power-like 
terms. It turns out that the longitudinal and transverse data cannot be explained simulta- 
neously by a renormalon contribution only; however, we will show that in the whole Q 2 range 
between 1 and 20 (GeV/c) 2 the longitudinal channel appears to be consistent with a pure 
Ji?-renormalon picture, adopting strengths which are not inconsistent with those expected 
in the naive non-abelianization (NNA) approximation |1| and agree well with the results of 
a recent analysis [|HJ of the CCFR data [|18| on xF^(x, Q 2 ) . Then, after including the 1/Q 2 
and 1/Q 4 J-R-renormalon terms as fixed by the longitudinal channel, the twist-4 and twist-6 
contributions arising from multiparton correlations are phenomenologically determined in 
the transverse channel and found to have non-negligible strengths with opposite signs. It is 
also shown that our determination of multiparton correlation effects in F^(X, Q 2 ) is only 
marginally affected by the specific value adopted for the strong coupling constant at the 
Z-boson mass, at least in the range between ~ 0.113 and ~ 0.118. Finally, an interesting 
outcome of our analysis is that the extracted twist- 2 contribution in the deuteron turns out 
to be compatible with the enhancement of the ci-quark parton distribution recently proposed 
in Ref. || . 



The paper is organised as follows. In the next Section the main features of the OPE, 
the NLO approximation of the leading twist and the usefulness of the Natchmann definition 
of the moments are briefly reminded. In Section 3 our procedure for the evaluation of the 
experimental longitudinal and transverse moments is described. Section 4 is devoted to a 
phenomenological analysis of the data at NLO, while the inclusion of the 7i?-renormalon 
uncertainties is presented in Section 5. The results of our final analysis of the transverse 
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data are collected in Section 6, and our main conclusions are summarised in Section 7. 



The Operator Product Expansion and the Leading 
Twist 



The complete Q 2 evolution of the structure functions can be obtained using the OPE [2p_ 
of the time-ordered product of the two currents entering the virtual-photon nucleon forward 
Compton scattering amplitude, viz. 



T[J(z) J(0)] 



E 



f*(-z 2 ) Z^Z^.-.Z^ 



O 



li±H2...ll n 



where 2 „ are symmetric traceless operators of dimension d% and twist r° = d^—n, with 
a labelling different operators of spin n. In Eq. ([I]) f^(—z 2 ) are coefficient functions, which 
are calculable in pQCD at short- distance. Since the imaginary part of the forward Compton 
scattering amplitude is simply the hadronic tensor containing the structure functions mea- 
sured in deep inelastic scattering (DIS) experiments, Eq. (fj) leads to the well-known twist 
expansion for the Cornwall-Norton {CN) moments of the transverse structure function, viz. 



M T n {Q 2 



x n ~ 2 F?(x,Q 2 ) 



£ E nT ^,Q 2 )O n M (gr 

T=2,even V^* 



(2) 



where /i is the renormalization scale, O nT (fi) are the (reduced) matrix elements of operators 
with definite spin n and twist r, containing the information about the non-perturbative 
structure of the target, and E nr (fi, Q 2 ) are dimensionless coefficient functions, which can be 
expressed perturbatively as a power series of the running coupling constant a s {Q 2 ). 

As it is well known, in case of the leading twist (r = 2) one ends up with a non-singlet 



quark operator 0^=2 with corresponding matrix elements O 



NS 



,NS 



and coefficients E, 



NS 
n2 ) 



and with singlet quark 0^ =2 and gluon 0^ =2 operators with corresponding matrix elements 
and coefficients E^ 2 (after their mixing under renormalization group equations); explicitly 
one has 



Ml{Q 2 ) = f^(Q 2 ) + higher twists 



(3) 
(4) 



with at NLO 



,NS 




l + a s (Q 2 )R^/4n 
l + a s (ii 2 )R^/4ir 

l + a s (Q 2 )R-/An 
1 + a s (// 2 )i?-/47r 

l + a s (Q 2 )R+/A7T 
l + a s {fi 2 )R+/A7T 



+ 



(5) 



(6) 
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where all the anomalous dimensions 7* 
Ref. 



NS, ±) and coefficients RL can be found in, 



e.g 



21 



for n < 10 and for a number of active flavours, Nf, equal to Nf = 3, 4 and 5. 
At NLO the running coupling constant a s (Q 2 ) is given by 



47T 



/5 ln(QVA 2 



1 - 



ft In ln(Q 2 /A 2 ; 



MS' 



Pi 



ln(Q 2 /A 2 



(7) 



MS' 



where Aj^g is the QCD scale in the MS scheme, /3q 
The coefficients or can be rewritten as follows 



11 - 2JV//3 and fa = 102 - 38JV//3. 



(1 - b n ) ^ 2 ) 
b n /^(/i 2 ) - c r 



+ c, 



n 

2\ 



G 



where fj%(Q ) is the CA moment of the gluon distribution function Q ) of order n, viz. 



fiffl = 1} dx a? 
large x, one has i> r , 



2 G(x,Q 2 ) 
^ 1 and c n 



It turns out that, since quark and gluons are decoupled at 
> already for n > 4, so that only a~ contributes to Eq. (|6|). 



Moreover, again for n > 4 one has 7 n ~ 7^ and R n (f) 



r>NS 
,S(n2- 



Bj], which implies that for 



n > 4 the evolution of the leading-twist singlet moments ) almost coincide with that 

of the non-singlet ones fi^ s (Q 2 ). Therefore, when n > 4, we will consider for the leading 
twist term the following NLO expression 



^n S {Q 2 )+^{Q 2 ) 



i>4 a { n 



a s (/i 2 



l + a s (Q 2 )Rn 3 /^ 
l + a s (/i 2 )^/4vr 



(9) 



with affl = fx^ s (fi 2 ) + (fJ 2 ). To sum up, the Q 2 evolution of the leading-twist transverse 
moments is completely determined by pQCD for each spin n and all the unknown twist-2 
parameters reduce to the three matrix elements a^ s and af in case of the second moment, 
and only to one matrix element, affl, for n > 4. 

The contribution of power-like corrections to the Q 2 dependence of the moments 
(0) is due to higher twists corresponding to r > 4. Already at r = 4 the set of basic 
operators is quite large and their number grows rapidly as r increases; moreover, each set of 
operators mix under renormalization group equations. The short- distance coefficients E nr 
can in principle be determined perturbatively, but the calculations are cumbersome and the 
large number of the involved matrix elements O nT (/j,) makes the resulting expression for the 
moments of limited use. Therefore, following Ref. |R|, in this work we will make use of 
effective anomalous dimensions for the higher-twist terms (see next Section). 

The OPE can be applied also to the longitudinal structure function Fff(x, Q 2 ), where 
power corrections are expected to be more important than in the transverse case, because 
F^(x, Q 2 ) is vanishing at LO. However, NLO effects generate a non- vanishing contribution; 
more precisely, since the longitudinal structure function is defined as 



F?(x,Q 2 ) = F 2 N (x } Q 2 ) 1 + 



4M¥ l 
Q 2 , 



R L / T [ 



x,Q 2 



1 + R% /T \ 



x,Q 2 



(10) 
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where R^/ T (x,Q 2 ) = a^(x,Q 2 )/a^(x,Q 2 ) is the L/T cross section ratio, one has 

M^(Q 2 ) = f 1 dx x n ~ 2 F?(x, Q 2 ) = Hn(Q 2 ) + ni gh er twists (11) 
Jo 

with the leading-twist contribution /i^(Q 2 ) given at NLO by 

Am = a f* ) n \ l {r im + ^h^ Q2) } + ^ m (12) 

where fiKQ 2 ) is the light-quark flavour contribution to the n-th moment of the transverse 
structure function, fi^(Q 2 ) is the charm quark contribution (starting when the invariant 
produced mass W is greater than twice the mass of the charm quark), and d = 12/9 and 
20/9 for Nf = 3 and 4, respectively. In what follows Eq. (|i~2"D will be evaluated using sets of 
parton distributions available in the literature (see Section 5). 

For massless targets only operators with spin n contribute to the n-th CN moments 
of both the longitudinal and transverse structure functions. When the target mass is non- 
vanishing, operators with different spins can contribute and consequently the higher-twist 
terms in the expansion of the CN moments M^(Q 2 ) and M^(Q 2 ) contain now also target- 



mass terms, which are of pure kinematical origin. It has been shown by Natchmann [22| that 
even when M/fl the moments can be redefined in such a way that only spin-n operators 
contribute in the n-th moment, namely 



M^(Q 2 ) = f dx 
Jo 



x 3 



F N (x O 2 ) + A -^F N (x 2 > + l)£/s-2(n + 2) 



(14) 



where r = y 1 + 4M 2 x 2 /Q 2 and £ = 2x/(l + r) is the Natchmann variable. Using the 
experimental FJ?(x, Q 2 ) and F^(x, Q 2 ) in the r.h.s. of Eqs. (|I~3-|i~4]), the target-mass correc- 
tions are cancelled out and therefore the twist expansions of the experimental Natchmann 
moments M^(Q 2 ) and M^(Q 2 ) contain only dynamical twists, viz. 

Mn(Q 2 ) — fJ'niQ 2 ) + dynamical higher twists (15) 
M^(Q 2 ) = ^{Q 2 ) + dynamical higher twists (16) 

where the leading-twist terms /j^(Q 2 ) and n^{Q 2 ) are given at NLO by Eqs. (§-|D and flljD , 
respectively. 



3 Data on Transverse and Longitudinal Moments 

For the evaluation of the Natchmann transverse Mj(Q 2 ) (Eq. ([HJ)) as well as longitudi- 
nal M^(Q 2 ) (Eq. (0)) moments systematic measurements of the experimental structure 
functions FJ?(x, Q 2 ) and F^(x, Q 2 ) are required in the whole x-range at fixed values of Q 2 . 



7 



However, such measurements are not always available and therefore we have adopted inter- 
polation formulae ("pseudo-data"), which fit the considerable amount of existing data on 
the proton and deuteron structure functions. 

In case of the transverse channel we have used: i) the Bodek's fit [§, [7| to the 
inclusive (e, e') SLAC data in the resonance region for values of the invariant produced 
mass W smaller than 2.5 GeV, and ii) the Tulay's fit || to the world data in the DIS 
region for W > 2.5 GeV. In order to have continuity at W = 2.5 GeV the Bodek's fit has 
been divided by a factor 1.03 (at any Q 2 ), which is inside the experimental errors. Moreover, 
we have also taken into account the wealth of SLAC proton data [§] beyond x ~ 0.7 up to 
x ~ 0.98 through a simple interpolation fit given in Ref. ||. All these fits cover the range of 
x which is crucial for the evaluation of the moments considered in this work; therefore, the 
uncertainties on the moments are related only to the accuracy of the interpolation formulae; 
the latter are simply given by a ±4% total (systematic + statistical) error reported for the 
SLAC data [Q, [J and by the upper and lower bounds of the Tulay's fit quoted explicitly in 
Ref. ||. Finally, since the whole set of DIS (unpolarised) SLAC and BCDMS data (on 
which also the Tulay's fit is based) is known to favour the value a; s (M§) ~ 0.113 (see Ref. 
|3), in what follows the value = 290 (240) MeV at N f = 3 (4) will be adopted for the 



calculation of the running coupling constant at NLO via Eq. 

As for the longitudinal channel, the structure function F^(x, Q 2 ) is reconstructed via 
Eq. (|ToD , using for the L/T ratio the phenomenological fit (and the interpolation uncertain- 
ties) provided in Ref. |TI[. In the nucleon resonance region the ratio R^ T (x,Q 2 ) has been 
taken from Ref. [24| for the lowest values of Q 2 (namely Q 2 = 1 4- 2 (GeV/c) 2 ), while it has 



been assumed to be equal to zero for Q 2 > 2 (GeV/c) 2 , which is consistent with the limited 



amount of available data |]25| . 

Since the OPE is totally inclusive, the contribution of the nucleon elastic channels 
has to be included in the calculation of the moments, viz. 



F»(x,Q'< 



f {inel \x,Q 2 )+F 2 N(el \x,Q 2 ) 
^ {inel \x,Q 2 )+F? (el \x,Q 2 ) 



(17) 



and consequently 

MJ(Q 2 ) 
M n L (Q 2 ) 

In case of the proton one has 



M T n {Q 2 ) 
M^(Q 2 ) 



inel 



+ 



inel 



+ 



M n(Q 2 ] 
M n(Q 2 ) 



F p{el) { 



x,Q 2 



x,Q' 



8{x 
S(x 



[GUQ 2 )} 2 + V [G p M (Q 2 )} 2 



1 + 7/ 



[GUQ 2 )? 

V 



b The impact of a different value of a s (M§ ), much closer to the updated PDG value a s {M\ 



(19) 

(20) 
0.119 ± 



0.002 



will be briefly addressed at the end of Section 5. 



s 



where G P E {Q 2 ) and G P M {Q 2 ) are the electric and magnetic (Sachs) proton form factors, respec- 
tively, and n = Q 2 /4M 2 . Thus, the contribution of the proton elastic channel is explicitly 
given by 



el 



proton 



M^{Q 2 ) 



el 



'proton 



n+1 



3 + 3(n + l)r* + n(n + 2)r 



*2 



l+r*J (n + 2)(n + 3) 

\GUQ 2 )f + v[G p M {Q 2 )] 2 



1 + n 

2 \ n+1 



l + r] V 1 + r* 
n + 1 
n + 3 



{[G p E (Q 2 )] 2 -[G p M (Q 2 )} 2 + 



(21) 



(n + 2)(l + r*) 



[G p E (Q 2 )} 2 + v[G p M (Q 2 )} 2 ' 
i + v 



(22) 



with r* = y/1 + AM 2 /Q 2 = y/l + 1/rj. 

In case of the deuteron the folding of the nucleon elastic channel with the momentum 
distribution of the nucleon in the deuteron gives rise to the so-called quasi-elastic peak. In the 
Q 2 -range of interest the existing data are only fragmentary and therefore we have computed 
F 2 D(ei) (x, Q 2 ) and pf (e °(x, Q 2 ) using the procedure of Ref. PB] , which includes the folding of 
the nucleon elastic peak with the nucleon momentum distribution in the deuteron as well as 
final state interaction effects; such procedure can be applied both at low (Q 2 < 1 (GeV/c) 2 ) 
and high values of Q 2 (see Ref. ||27||). The results of the calculations, performed using 
the deuteron wave function corresponding to the Paris nucleon- nucleon potential f2~8| , are 
compared in Fig. 1 with the SLAC data of Ref. [|| at Q 2 = 1.75 and 3.25 (GeV/c) 2 . It 
can be seen that the agreement with the data is quite good for the transverse part and still 
acceptable for the less accurate longitudinal response. 

In Figs. 2 and 3 the Q 2 behaviours of the inelastic, elastic and total transverse 
moments are reported for the proton and the deuteron, respectively, in the whole Q 2 -range 
of interest in this work, i.e. 1 < Q 2 {GeV/ c) 2 < 20. Note that in case of the deuteron our data 
refer to the moments of the structure function per nucleon. It can clearly be seen that the 
inelastic and elastic contributions exhibit quite different Q 2 -behaviours and, in particular, 
the inelastic part turns out to be dominant for Q 2 > 1 (GeV/c) 2 in the second moment and 
only for Q 2 > n (GeV/c) 2 in higher order moments. The Q 2 -behaviour of the longitudinal 
moments is illustrated in Figs. 4 an 5 in case of the proton and the deuteron, respectively. 
It can be seen that the elastic contribution drops quite fast and it is relevant only at the 
lowest values of Q 2 (Q 2 ~ few (GeV/c) 2 ). 

We point out that our present determination of the longitudinal and transverse ex- 
perimental moments is clearly limited by the use of phenomenological fits of existing data 
(i.e. "pseudo-data"), which are required in order to interpolate the structure functions in 
the whole x-range for fixed values of Q 2 , as well as by the existing large uncertainties in 
the determination of the L/T ratio R^ T (x,Q 2 ). Thus, both transverse data with better 
quality at x > 0.6 and Q 2 < 10 (GeV/c) 2 and more precise, systematic determinations of 
the L/T cross section ratio are still required to improve our experimental knowledge of the 
Q 2 -behaviour of the low-order moments of the nucleon structure functions. 
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4 Analysis of the Transverse Data at NLO 



In this Section we present our analysis of the transverse pseudo-data of Figs. 2 and 3, 
adopting for the leading twist the NLO expressions and for the power corrections 

a purely phenomenological ansatz. Indeed, as already pointed out, several higher-twist 
operators exist and mix under the renormalizat ion-group equations; such a mixing is rather 



involved (see, e.g., Ref. |30[ in case of twist-4 operators) and in particular the number of 
mixing operators increases with the spin n. Since complete calculations of the higher-twist 
anomalous dimensions are not yet available, we use the phenomenological ansatz already 
adopted in Ref. which in case of the second moment Mj =2 (Q 2 ) leads to the expansion 



M^Q 2 ) = ^ s (Q 2 ) + ^(Q 2 ) + a? 



a s {Q 2 
a s (fi 2 



(4) 

12 u 2 

w (23) 



where the logarithmic pQCD evolution of the twist-4 contribution is accounted for by the 
term [a s (Q 2 )/a s (/i 2 )] 7 2 with an effective anomalous dimension and the parameter a\ 
represents the overall strength of the twist-4 term at Q 2 = fi 2 . For the QCD renormalization 
scale /i we adopt hereafter the value /i = 1 GeV/ c. We have explicitly checked that a different 
choice for fi does not modify any conclusion of this work, because (e.g. in Eq. (|2"3"D) the 
values of the parameter a\ corresponding to two different choices [i and // turn out to be 
related by the logarithmic pQCD evolution, i.e. af(//) = a^/x) [a s {ji' 2 ) / a s {fi 2 )] 12 fi 2 / n' 2 . 

The unknown parameters appearing in Eq. ( p3[) (i.e., the three twist-2 parameters 
a^ s = /4 fS {Q 2 = A* ), af = fiHQ 2 = fJ 2 ) and aff = ^ (Q 2 = fj 2 ) and the two twist-4 
parameters cffi and have been determined by fitting our data for the proton and the 
deuteron (see Fig. 2(a) and 3(a), respectively), adopting the least-x 2 procedure in the Q 2 - 
range between 1 and 20 (GeV/c) 2 . It turned out that the singlet and gluon twist-2 parameters 
can be directly determined form the data only in case of the deuteron, but not in case of 
the proton. Thus, in fitting the proton data we have kept fixed the parameters af and 
at the values obtained for the deuteron. Our results are shown in Figs. 6(a) and 7(a) in 
the whole Q 2 -range of the analysis, while the values obtained for the various parameters are 
reported in Table 1, together with the uncertainties of the fitting procedure corresponding 
to one- unit increment of the x 2 /N variable (where N is the number of degrees of freedom). 
It can clearly be seen that: 

• the twist-4 contribution to second transverse moment M^iQ 2 ) is quite small in the 
proton and is almost vanishing in the deuteron in the whole Q 2 -range of our analysis; 
the latter result suggests that the twist-4 effect in the neutron comes with a sign 
opposite to that in the proton at variance with the expectations from the bag model 

• the twist-4 contribution in the proton at Q 2 = jj, 2 = 1 (GeV/c) 2 (i.e., a$ = 0.012 ± 
0.010) turns out to be significantly smaller than the result t4 = 0.030 ± 0.003 quoted 



in Ref. ||13|| , where however the twist-2 term was not simultaneously fitted to the 



data, but instead it was calculated using the parton densities of Ref. p2[ evolved at 
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NLO. Therefore, we have repeated our analysis by fixing the twist-2 term at the GRV 
prediction |H| obtaining aj> — 0.02 ± 0.01, which is still lower but not inconsistent 
within the errors with the result of Ref. ||13|| . We point out that a small variation of 
the twist-2 term can affect significantly the strength of the small twist-4 term; that is 
why our uncertainty from a simultaneous a^ 5 and fit appears to be quite larger 
than the one found in Ref. ||13|| : 

• as a cross-check, we have also fitted separately the non-singlet parameter a^ s to the 
"(p — n)/2" data, defined as the difference between the proton and deuteron data, 
obtaining a% s = 0.029 ± 0.009, = 0.012 ± 0.010 and = 1 ± 1. Combining 
these results with those found in the deuteron, we expect to get in case of the proton 
a% s = 0.096 ± 0.017, ajj 4) = 0.012 ± 0.010 and 7$ 4) = 1 ± 1, which indeed are in nice 
agreement with the proton fit results of Table 1. 

As described in Section 2, in case of the transverse moments M^(Q 2 ) with n > 4 the 
evolution of the twist-2 contribution can be simplified and assumed to be a pure non-singlet 
one; therefore, we have considered the following twist expansion 



Mr >4 (Q 2 )=/iI(Q 2 )+ai 4 ) 



a s (Q 2 



(4) 

T™ 2 

Q 2+ n 



a s (Q 2 
a s (fi 2 



$ (24) 



where the leading twist term /i^(Q 2 ) is given at NLO by Eq. (|9]) and fi = 1 GeV/c. All 
the five unknown parameters (i.e., and 7^ 6 * ) ) have been determined from 

the data through the least-x 2 procedure for each value of n. The results are shown in Figs. 
6 and 7, while the values of the parameters are reported in Tables 2 and 3 in case of the 
proton and deuteron, respectively. Our main results can be summarised as follows: 

• our twist-2 term, extracted from the proton data together with the twist-4 and twist-6 
contributions, differs only slightly from the predictions obtained using the set of parton 
distributions of Ref. |32| evolved at NLO (see dashed lines in Fig. 6). We have checked 
that, by repeating our analysis with the twist-2 term fixed at the GRV prediction, the 
values of all the higher-twist parameters appearing in Eq. (^4|) change for n > 2 only 
within the errors reported in Tables 2 and 3; 

• the twist-4 and twist-6 contributions result to have opposite signs and, moreover, at 
Q 2 ~ 1 (GeV/c) 2 they are approximately of the same order of magnitude. The nega- 
tive sign of the twist-6 is clearly due to the fact that a fit with a twist-4 term alone 
overestimates the low Q 2 data (Q 2 ~ few (GeV/c) 2 ). We stress that the opposite 
signs found for the twist-4 and twist-6 terms make the total higher-twist contribution 
smaller than its individual terms; in particular, at large Q 2 the sum of the twist-4 and 
twist-6 contributions turns out to be a small fraction of the twist-2 term (< 10% for 
Q 2 >n(GeV/c) 2 ). 

• the values of the effective anomalous dimensions 7^ and 7^ for n = 4, 6, 8 result 
to be around 4.0 and 2.5, respectively, i.e. significantly larger than the values of the 
corresponding twist-2 anomalous dimensions (7^ s — 0.8 -r- 1.2 for n = 4, 6, 8 ||21|| ); 
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• the uncertainties on the different twist contributions due to the parameter fitting pro- 
cedure are always within ±15% (see Fig. 7 in case of the deuteron); 

• the twist expansions and (pip appear to work quite well for values of Q 2 down to 
~ 1 {GeVjcf. 

An interesting feature of our analysis of the transverse moments is that the leading- 
twist contribution is extracted from the analysis and not fixed by calculations based on 
a particular set of parton distributions. The comparison of our extracted twist-2 term 



with the predictions based on the GRV parton distributions [^| is shown in Fig. 8 for 
Q 2 > 5 (GeV/c) 2 . It can clearly be seen that our results and the GRV predictions agree 
quite well in case of the proton, whereas they differ significantly in case of the deuteron 
for n > 2. The inclusion of a new empirical determination of the nuclear effects in the 
deuteron, obtained in Ref. from the nuclear dependence of the SLAC data, increases 



only a little bit the disagreement for n > 2 (see dashed lines in Fig. 8). Since moments with 
n > 2 are mostly sensitive to the high- a; behaviour of the structure function, the question 
arises whether our extracted twist-2 moments can be explained by an enhancement of the 
(i-quark parton distribution, like the one advocated in Ref. [|TJ|, which explicitly reads as 
d(x) = d(x) + 0.1x(l + x)u(x). The modified GRV predictions, including also the empirical 
nuclear corrections for the deuteron, are shown by the solid lines in Fig. 8 and agree quite 
well with our results both for the proton and the deuteron. Therefore, our extracted twist- 
2 moments are clearly compatible with the hypothesis of an enhancement of the <i-quark 
distribution at large x. 

So far the power corrections appearing in Eqs. (^3|-|4]) have been derived assuming 



for the leading twist the NLO in the strong coupling constant. Since our main aim is to try 
to estimate the target-dependent power corrections generated by multiparton correlations, 
it is necessary to estimate the possible effects of higher orders of the perturbative series, 
which defines the twist-2 coefficient functions E^^.Q 2 ) appearing in Eq. (0). However, 
it is well known that such series is asymptotic at best and affected by the so-called IR- 
renormalon ambiguities P3J. More precisely, the only way to interpret consistently the 



large order behaviour of the pertubartion theory leads unavoidably to ambiguities which 
have the general form of power-suppressed terms. This happens because certain classes 
of high-order radiative corrections to the twist-2 (the so-called fermion bubble chains) are 
sensitive to large distances (i.e., to the non-perturbative domain). Nevertheless it should 
be stressed that the Ji?-renormalon contribution to the twist-2 has to cancel against the 
ultraviolet quadratic divergencies of twist-4 operators (see, e.g., Refs. [|4], |35|, [TBI). This 



means that each twist in the expansion (0) is not defined unambigously, while the entire 
sum (i.e., the complete calculation) is free from ambiguities. Therefore, while the data 
cannot be affected by 7i?-renormalon uncertainties, our fitting procedure, based on the 
separation among various twists and on the theoretical perturbative treatment of the leading 
twist only, can suffer from 7i?-renormalon ambiguities. The latter are target-independent 
quantities, being of pure perturbative nature, while genuine higher-twist effects are related 
to multiparton correlations in the target (i.e., target-dependent quantities). 

The /.R-renormalon picture has been applied to the phenomenology of deep inelastic 
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lepton-hadron scattering as a guide to estimate the x dependence of power corrections 
|15| . Such an estimate has been found to be a good guess in case of the proton and the 
deuteron structure functions and this fact may be understandable in terms of the notion of a 
universal Ji?-finite effective strong coupling constant |36| or in terms of the hypothesis of the 
dominance of the quadratically divergent parts of the matrix elements of twist-4 operators 



37 



Before closing this Section, we point out that the Ji?-renormalon contribution behaves 
as a power-like term, but it is characterised by twist-2 anomalous dimensions. Thus, our ob- 
servation that for n = 4, 6, 8 the effective anomalous dimensions 7^ 4 - ) and 7^ 6 - ) extracted from 
our NLO fit to the transverse data are significantly different from the corresponding twist-2 
anomalous dimensions, might be an indication of the presence of multiparton correlation ef- 
fects (at least) in the transverse channel. In order to try to disentangle the latter from large 
order perturbative effects, we will consider explicitly in the next Section the power-like terms 
associated to the IR renormalons as the general uncertainty in the perturbative prediction 
of the twist-2 |3F 



5 IR Renormalons and the Analysis of the Longitudi- 
nal Channel 

Within the naive non-abelianization (NNA) approximation the contribution of the renor- 
malon chains (i.e. the sum of vacuum polarisation insertions on the gluon line at one-loop 
level) to the non-singlet parts of the nucleon structure functions F^(x,Q 2 ) and F 2 N (x,Q 2 ) 
is given by |L4 



F^Q 2 ) 
Fi R (x,Q 2 ) 



f'dzF^^Q 2 

J x Z 



dzF 2 LT (^,Q 2 



a irD2(z) 



Q 2 
,C 2 (z) 



A 



Q 2 



Q 4 



(25) 
(26) 



where the constants A 2 and A[ R are related to the log-moments of an Ji?-finite effective 
strong coupling constant j36|, F[ T and F 2 LT are the leading- twist structure functions and the 
coefficient functions C%u\ and D^u) are given explicitly in Ref . |TJ| . Thus, the 7i?-renormalon 
contribution to the transverse and longitudinal moments is explicitly given by 



I Q 2 



A IR ~ 



(27) 



A^ 

Q 4 



Qtt 



(A? 
\ Q 2 

Sa s (Q 2 )D 4 (n) 
67r n + 1 



8a s (Q 2 )D 2 (n) An 



An 



n+1 n+2 
n + 1 



+ 



n + 3 



(21 
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where in the last equation we have used the non-singlet part of the NLO relation (|i~2|) . The 
coefficients C 2 ( 4 ) (n) and L>2{4) (n) read as follows (14 



C 2 {n) = - n -8 + - + ^— + -^- + 4S n (29) 

n n + 1 n + 2 

<5 4 (n) = ^ 2 -|n + 16-l-^--^--45 n 
2 2 nn + 1 n + 3 

L> 2 (n) = -n-4+- + —^- + — ^- + 4S n (30) 

n n + 1 n + 2 

~ , , 1 , 5 4 4 12 

D 4 (n) = -n 2 + -n + 8 — - — — - AS n 

2 2 n n + 1 n + 3 

(31) 

where S n = J2l=i As already mentioned, Eqs. (|27| - p8D describe the contributions of /i? 

renormalons to the non-singlet structure functions, while the more involved case of the singlet 



parts of the DIS structure functions has been recently investigated in Ref. [p9|| . There it 
has been shown that the difference between the 7i?-renormalon contributions to the singlet 
and non-singlet moments is not relevant for n > 4 thanks to the quark-gluon decoupling at 
large x. Therefore, for n > 4 it suffices to consider Eqs. (127H28D after substituting n^ s {Q 2 ) 
with nl{Q 2 ), given at NLO by Eq. (§. 



An interesting feature of the /.R-renormalon terms ( P7| - Pq ) is that they are mainly 
governed by the values of only two (unknown) parameters, and A{ R , which appear 
simultaneously both in the longitudinal and transverse channels at any value of n. The signs 
of A 1 ^- and A™ are not theoretically known, because they depend upon the prescription 
used to circumvent the renormalon singularities of the Borel integrals, while within the 



NNA approximation their absolute values may be provided by |L5) 

f^A M 5/3 ir 10/3 (32) 

1 2 1 33-2^/ ' 1 4 1 33-2^/ K ' 

with Cp = 4/3. Using for A-^ the same values adopted for the NLO calculation of a s {Q 2 ) 
(see previous Section), one gets that l^^l varies from 0.10 to 0.13 GeV 2 , while l-A^I ranges 
from 0.015 to 0.030 GeV 4 for N f = 3,4. 

First of all, we have checked whether the power corrections to the NLO twist-2 
contribution in the transverse channel can be explained by pure 7i?-renormalon terms in the 
whole range 1 < Q 2 (GeV/c) 2 < 20. It turns out that: i) the quality of the resulting 
fit is not as good as the one obtained via the expansion (^) and the obtained minima of 
the x 2 /N variable can be much larger than 1; ii) the extracted values of l^^l an d 1^4^! 
result to be much greater than the NNA expectations (|32"D and to depend strongly upon 
the inclusion of the data at low Q 2 (Q 2 ~ few (GeV/c) 2 ). Moreover, the Ji?-renormalon 
contribution to the longitudinal moments (Eq. (^)), calculated using the values of A 2 R and 
A 1 / 1 determined from the analysis of the transverse channel, leads to a large overestimation of 
the longitudinal data, as already noted in Ref. |]I5| . These results may be viewed as an effect 
of the presence of higher-twist terms generated by multiparton correlations in the transverse 
data for Q 2 > 1 (GeV/c) 2 . To make this statement more quantitative, we start with the 
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analysis of the longitudinal data adopting for the power corrections a pure Ji?-renormalon 
picture; namely, for n > A we have used the following expansion 



M L n {Q 2 



^(Q 2 )+^ IR) (Q 2 ) 



>n>4 



^{Q 2 )+^{Q 2 



A? 
Q 2 



8a s (Q 2 ) D 2 (n) An 



6/T n + 1 n + 2 



+ 



A? 

Q 4 



8a s (Q 2 )D 4 (n) 
67r n + 1 



An 



n + 1 



n 



(33) 



where n^{Q 2 ) is given by Eq. ( |i~2"D and calculated using the G-RV parton distributions [[32 
while fi^(Q 2 ) (see Eq. (^)) is taken from the NLO analysis of the transverse data made in 
the previous Section (see Tables 2 and 3 for the values of the twist-2 parameters affl). 

For each value of n > A we have determined the values of A^ R and A 1R from the 
least-% 2 fit to the longitudinal data in the whole range 1 < Q 2 (GeV/c) 2 < 20; our results 
are reported in Tables 4-5 and Figs. 9-10 in case of proton and deuteron targets, respectively. 
It can clearly be seen that: i) the extracted values of A% R are almost independent of n, 
which means that the n-dependence (i.e., the shape in x) of the 1/Q 2 power correction is 
nicely predicted by the NNA approximation; ii) the values obtained for \A\ R \ are only 
slightly larger than the NNA expectation (0); iii) the determination of A{ R is almost 
compatible with zero and affected by large uncertainties, since the 1/Q 4 power corrections 
turn out to be quite small in the longitudinal channel; nevertheless; the extracted values 
are not completely inconsistent with the NNA predictions (j3~2D ; iv) the power corrections 
appear to be approximately the same in the proton and deuteron longitudinal channels (as 
it can be expected from a pure 7i?-renormalon phenomenology). 

As for the second moment M^iQ 2 ) we have simplified our analysis by taking only the 
non-singlet 7i?-renormalon contribution (|28|) , i.e. by totally neglecting its singlet part. This 
is an approximation, but the resulting fit to the data on Mf(Q 2 ) turns out to be quite good 
as it can be seen from Figs. 9(a) and 10(a), yielding values of A^ R and A^ R only slightly 
different from the ones previously determined by the analysis of the moments with n > A 
(see Tables 4 and 5). 

To sum up, in case of both proton and deuteron targets a pure 7i?-renormalon de- 
scription of power corrections works quite nicely in the longitudinal channel starting already 
at Q 2 ~ 1 (GeV/c) 2 . Averaging the results of Tables 4 and 5 for n > A only, our determina- 
tion of the 7i?-renormalon strength parameters results to be: A 2 — —0.132 ± 0.015 GeV 2 
and A\ R ~ 0.009 ±0.003 GeV A , which, we stress, are not inconsistent with the NNA expec- 
tations (|3"2"[). Moreover, the value found for A 1 ^ 1 nicely agrees with the corresponding findings 
of Ref. 0, recently obtained from a NLO analysis of the CCFR data pf on xF i f ( x , Q 2 )- 

Before closing this section, we stress again that the Ji?-renormalon power-like terms 
should be regarded as the general uncertainty in the perturbative calculation of the twist-2 
term. Thus, it is worth recalling that our estimates of the Ji?-renormalon parameters have 
been obtained by taking the perturbative calculation of the twist-2 term at NLO. We expect 
that our determination of the Ji?-renormalon parameters holds only at NLO and would vary 
if higher orders of the perturbation theory were included. As a matter of fact, a significative 
reduction of the 7i?-renormalon terms seems to be suggested by the recent NNLO results 
quoted in Refs. |16[ and flTT|. 
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6 Final Analysis of the Transverse Data 



After having determined the strengths and signs of the twist-4 and twist-6 Ji?-renormalon 
contributions from the analysis of the longitudinal channel, we can now proceed to the final 
analysis of the transverse data, adopting the following twist expansion 



M T n {Q 2 )=f n {Q 2 )+^ IR \Q 2 ) + a 



(4) 



a s (Q 2 



7n 



^- + a (6) 
Q2 +a n 



a s (/i 2 



W 6) 

in 



Q 4 



(34) 



where now the higher-twist terms involving the parameters affl , 7^ , and 7^ should be 
related to (target-dependent) multiparton correlation effects. Collecting Eqs. @ and 
one has for n > 4 



M n T >4 (Q 2 



,(2) 



a s (n 2 ) 



In 



In 



l + a s (Q 2 )R^ s /4n 


{" 


A? 




Q 2 


^- + a (6) 


[« S (Q 2 )1 

_a s (/i 2 )_ 


In 





C 2 (n) 



A IR ~ 1 



Q 



(35) 



The Ji?-renormalon parameters A 2 R and A4 have been kept fixed at the values A 2 R = 
—0.132 GeV 2 and A IR = 0.009 GeV 4 found in the previous Section, while the values of 
the five parameters a%\ affl, 7^', a^' an d have been determined through the least- 
X 2 procedure and reported in Tables 6 and 7 in case of proton and deuteron, respectively. 
Comparing with the results obtained without the /.R-renormalon contributions (see Tables 
2 and 3), it can be clearly seen that the values of the twist- 2 parameters affl are almost 
unchanged, while the values of the higher-twist parameters af£\ 7^, and 7^ vary only 
within the uncertainties of the fitting procedure. Note that with the inclusion of the IR- 
renormalon contribution the sum a$ + affl is closer to zero, which implies that at Q 2 ~ 
/i 2 = 1 (GeV/c) 2 the twist-4 and twist-6 terms generated by multiparton correlation almost 
totally compensate each other. Such an effect is clearly illustrated in Figs. 11 and 12, 
where the contributions of the twist- 2 at NLO, of the IR- renormalon and of the multiparton 
correlations are separately reported. From Figs. 11 and 12 it can be also seen that, for n > 4, 
the IR- renormalon contribution increases significantly around Q 2 ~ 1 (GeV/c) 2 and could 
become of the same order of magnitude of the twist- 2 term at NLO in case of higher order 
moments. The effects from multiparton correlations appear to exceed the 7i?-renormalon 
term only for Q 2 > 2 (GeV/c) 2 (at n > 4). 

As for the second moment M 2 {Q 2 ), following our previous analyses, we apply the IR- 
renormalon correction only to the non-singlet twist-2 term; moreover, the twist- 2 parameters 
af and a 2 are kept fixed at the values given in Table 1 for the deuteron and only the twist-4 

(4) 

term a$ [a s (<5 2 )/a s (/i 2 )] 7n {is 2 /Q 2 ) is explicitly considered in the analysis. The resulting 
value of the twist-2 parameter a 2 s is 0.096±0.006 (0.066±0.0044) for the proton (deuteron), 
which coincides within the uncertainties with the one given in Table 1. The twist-4 parameter 



a 2 ^ 
and 



turns out to be almost compatible with zero, namely a^f 1 = 0.01 ± 0.01 for the proton 
\ a ^\ %j 10~ 3 f° r the deuteron. Considering also the results of the previous Section on 
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the longitudinal channel, our analyses indicate that the smallness of multiparton correlation 
effects on the second moments (both transverse and longitudinal ones) is consistent with the 
Q 2 -behaviour of the data starting already at Q 2 

H 



Basing on naive counting arguments (see, 



: 1 {GeV/cf 
e.g., Ref. 



one can argue that the 



twist expansion for the transverse moments at Q ~ 11 can be approximately rewritten as 



<G" 2 ) - Ag> 



1 + n 



71' 



"r( 6 )' 

n 



(36) 



where is the twist-2 contribution and T^ 4 ^ (T^) represents the mass scale of the twist-4 
(twist-6) term, expected to be approximately independent of n for n > 4. (Note that in Eq. 
(Jj6|) we have already taken into account the opposite signs of the twist-4 and twist-6 terms 
as resulting from our analyses). Thus, one gets 



F (4) 



\ 



(4) 

ah 



nA 



(2) 



F (6) 



,.(6)1 



n 2 A 



(2) 



1/4 



(37) 



Our results for F^' and r^ 6 ) at n = 4, 6, 8, obtained taking affl and from Tables 6 and 7 
and using for Affl the twist-2 term at NLO (i.e., a^) plus the whole 7i?-renormalon contri- 
bution (as determined from our fitting procedure and evaluated at Q 2 = fi 2 = 1 (GeV/c) 2 ), 
are collected in Fig. 13. It can clearly be seen that the mass scales of the twist-4 and 
twist-6 terms are indeed approximately independent of n, viz. 

r (4) ^ r (4) ^ q 7 q GeV and 

r {6) ^ r (6) ^ Q 55 GeV The yalue obtained for p(4) 

is significantly higher than the naive ex- 
pectation r^ 4 ^ ~ J < k 2 ± > ~ 0.3 GeV |40], p]] , but not very far from the results of Ref. Ul3"| . 



Without including the Ji?-renormalon contribution in A^ (i.e., taking only Affi 



,.(2) 



the 

values of I^ 4 ) and r( 6 ) would increase by ~ 20% and ~ 10%, respectively (cf. also Ref. |42| ). 

Before closing the Section, we want to address briefly the impact that the specific value 
adopted for the strong coupling constant at the Z-boson mass can have on our determination 
of multiparton correlation effects in the transverse channel. As already mentioned, existing 
analyses of the whole set of DIS (unpolarised) world data favour the value a s (M|) ~ 0.113 
Ref. [12]), which is however well below the updated PDG value a s (M|) = 



(see, e.g., 

0.119 ± 0.002 p3]| . Moreover, in Ref. [0 it has been argued that an increase of a s (M|) up 
to ~ 0.120 can give rise to a significative decrease of the relevance of the higher-twists effects 
in the DIS data (up to a reduction by a factor ~ 2). Therefore, we have repeated our analyses 
of longitudinal and transverse pseudo-data adopting the value ct s {M\ ) = 0.118, where a set 
of parton distributions is available from Ref. j|3|| (we use the parton distributions only for 
the calculation of li^(Q 2 ) at NLO). All the results obtained at the higher value of a s (M|) 
have the same quality as those presented at a s (M|) = 0.113 with a slight, but systematic 
increase of the minima of the x 2 /N variable. The 7i?-renormalons parameters A^ R and A\ R 
are determined again by fitting the longitudinal data in the Q 2 -range from 1 to 20 (GeV/c) 2 , 
and their values are now given by A% — —0.103 ±0.017 and A\ R ~ 0.005 ±0.004, which are 
compatible within the quoted uncertainties with the corresponding results of the high-Q 2 
analysis of Ref. [fTEfl. Thus, by construction, the sum of the twist-2 term at NLO and the 
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IR- renormalon contribution is almost independent of the specific value of a g (M§) in the 
longitudinal channel. However, the same happens in the transverse channel as it is clearly 
illustrated in Fig. 14, where the results obtained at a s (M|) = 0.113 and a; s (M§) = 0.118 
are compared in case of the transverse moments with n > 4 and found to differ only by 
less than ~ 5%. Therefore, though the NLO twist-2 terms as well as the 7i?-renormalon 
contributions are separately sensitive to the specific value of a s (M|), their sum turns out 
to be quite independent of q; s (M|). This means that our determination of the multiparton 
correlation effects is only marginally affected by the specific value adopted for a s (M|), at 
least in the range of values from ~ 0.113 to ~ 0.118. 



7 Conclusions 

We have analysed the power corrections to the Q 2 behaviour of the low-order moments 
of both the longitudinal and transverse structure functions of proton and deuteron using 
available phenomenological fits of existing data in the Q 2 range between 1 and 20 (GeV/c) 2 . 
The SLAC proton data of Ref. ||, which cover the region beyond x ~ 0.7 up to x ~ 0.98, 
as well as existing data in the nucleon-resonance regions have been included in the analysis 
with the aim of determining the effects of target-dependent higher- twists (i.e., multiparton 
correlations) . 

The Natchmann definition of the moments has been adopted for disentangling prop- 
erly kinematical target-mass and dynamical higher-twist effects in the data. The leading 
twist has been treated at the NLO in the strong coupling constant and, as far as the 
transverse channel is concerned, the twist-2 has been extracted simultaneously with the 
higher-twist terms. The effects of higher orders of the perturbative series have been esti- 
mated adopting the infrared renormalon model of Ref. fll4] |, containing both 1/Q 2 and 1/Q 4 



power-like terms. It has been shown that the longitudinal and transverse data cannot be 
explained simultaneously by the renormalon contribution only; however, in the whole Q 2 
range between 1 and 20 (GeV/c) 2 the longitudinal channel appears to be consistent with a 
pure Ji?-renormalon picture, adopting strengths not inconsistent with those expected in the 
naive non-abelianization approximation and in nice agreement with the results of a recent 
analysis |T7| of the CCFR data [|1^] on xF^(x, Q 2 ). Then, after including the 1/Q 2 and 1/Q 4 



Ji?-renormalon terms as fixed by the longitudinal channel, the contributions of multiparton 
correlations to both the twist-4 and twist-6 terms have been phenomenologically determined 
in the transverse channel and found to have non-negligible strengths with opposite signs. It 
has been also checked that our determination of the multiparton correlation effects is only 
marginally affected by the specific value adopted for a s (M|) (at least in the range from 
~ 0.113 to ~ 0.118). 

An interesting outcome of our analysis is that the extracted twist-2 contribution 
in the deuteron turns out to be compatible with the enhancement of the ci-quark parton 
distribution recently proposed in Ref. [|T^ . 

Let us stress that our analysis is presently limited by: i) the use of phenomenolog- 
ical fits of existing data (i.e., pseudo-data), which are required in order to interpolate the 
structure functions in the whole arrange for fixed values of Q 2 , and ii) by the existing large 
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uncertainties in the determination of the L/T ratio R^ T (x, Q 2 ). Therefore, both transverse 
data with better quality at x > 0.6 and Q 2 < 10 (GeV/c) 2 and more precise and systematic 
determinations of the L/T cross section ratio, which may be collected at planned facilities 
like, e.g., J LAB @ 12 GeV, could help to improve our understanding of the non-perturbative 
structure of the nucleon. Finally, we want to point out that, since in inclusive data multi- 
parton correlations appear to generate power-like terms with opposite signs, seminclusive or 
exclusive experiments might offer the possibility to achieve a better sensitivity to individual 
non-perturbative power corrections. 
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TABLES 



Table 1. Values of the twist-2 parameters a^ s = /x^ 5 (/x 2 ), af = ^f (yU 2 ) and af = /z[?(/z 2 ), 
and of the twist-4 parameters a^ and 7^ 4 \ obtained by fitting with Eq. (^) the pseudo-data 
for the proton and the deuteron (see Figs. 2(a) and 3(a), respectively). In case of the proton 
the parameters af and aSf have been kept fixed at the values obtained for the deuteron. The 
last row reports the value of the x 2 °f the fit divided by the number of degrees of freedom. 



Table 1 



parm 


deuteron 


proton 




0.067 ±0.014 


0.092 ±0.007 


af 


0.112 ±0.015 


0.112 


aSf 


0.066 ± 0.040 


0.066 


a (4) 
a 2 


< 10" 4 


0.012 ±0.010 


(4) 
72 




1 ± 1 


X 2 /N 


0.11 


0.10 



Table 2. Values of the twist-2 parameter and of the higher- twist parameters 7^, 
a$ and 7^, obtained by fitting with Eq. ( fM|) the proton pseudo-data of Fig. 2 for n = 4, 6 
and 8. The last row reports the value of the \ 2 °f the fit divided by the number of degrees 
of freedom. 



Table 2 



parm 


proton 


Mi 


Mi 


Mi 




0.028 ±0.001 


0.0078 ±0.0003 


0.0030 ±0.0001 




0.071 ±0.009 


0.046 ±0.003 


0.029 ±0.003 


in 


3.8 ±0.4 


3.7 ±0.2 


3.6 ±0.2 


a n 


-0.056 ±0.008 


-0.039 ±0.005 


-0.026 ±0.002 


~(6) 
in 


2.5 ±0.5 


2.3 ±0.2 


2.1 ±0.4 


X 2 /N 


0.05 


0.27 


0.65 
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Table 3. The same as in Table 2, but for the deuteron. 



Table 3 



parm 


deuteron 


M[ 


Mi 


Mi 




0.022 ±0.001 


0.0060 ±0.0002 


0.0024 ±0.0001 


u n 


0.055 ±0.006 


0.036 ±0.004 


0.022 ±0.002 


In 


3.9 ±0.2 


4.0 ±0.2 


3.8 ±0.3 


u n 


-0.046 ± 0.005 


-0.032 ±0.004 


-0.020 ±0.002 


~(6) 
In 


2.5 ±0.2 


2.4 ±0.2 


2.2 ±0.5 


X 2 /N 


0.10 


0.37 


1.0 



Table 4. Values of the Ji?-renormalon parameters A^ R and A^ R obtained using Eq. to 
fit the longitudinal proton data of Fig. 4. The twist-2 contribution /i^(Q 2 ) is calculated at 
NLO via Eq. (|i~2"D adopting the GRV set of parton distributions f32|, while fj,^(Q 2 ) (see Eq. 
(§)) is obtained using the values of aff reported in Table 2. The last row reports the value 
of the x 2 °f the fit divided by the number of degrees of freedom. 



Table 4 



parm 


proton 


M£ 


M£ 




K 


Ai R 


-0.16 ±0.05 


-0.12 ±0.04 


-0.11 ±0.04 


-0.13 ±0.04 


A 1R 


0.06 ±0.03 


0.01 ±0.01 


-0.001 ±0.001 


0.0004 ±0.0004 


X 2 /N 


0.24 


1.5 


1.1 


0.42 
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Table 5. The same as in Table 4, but for the deuteron. 



Table 5 



parm 


deuteron 




M£ 


M£ 


M- 


A{* 


-0.18 ±0.06 


-0.14 ±0.04 


-0.14 ±0.04 


-0.15 ±0.01 


A* 


0.06 ±0.03 


0.025 ±0.013 


0.012 ±0.008 


0.009 ±0.003 


X 2 /N 


0.09 


0.26 


0.36 


0.18 



Table 6. Values of the twist- 2 parameter affi and of the higher-twist parameters af£\ 
7^ 4 - ) , and ^\ obtained by fitting with Eq. ( j35|) the proton pseudo-data of Fig. 2 for 
n = 4, 6 and 8. The Ji?-renormalon twist-4 and twist-6 contributions are obtained using 
A[ R = -0.132 GeV 2 and A{ R = 0.009 GeV 4 (see text). The last row reports the value of 
the x 2 °f the fit divided by the number of degrees of freedom. 



Table 6 



parm 


proton 


Mi 


Mi 


Mi 




0.028 ±0.001 


0.0079 ±0.0003 


0.0031 ±0.0001 


a n 


0.061 ±0.008 


0.045 ±0.002 


0.028 ±0.003 


~(4) 

in 


4.1 ±0.2 


4.1 ±0.2 


4.1 ±0.5 




-0.051 ±0.008 


-0.041 ±0.002 


-0.027 ±0.003 


in 


2.7 ±0.2 


2.6 ±0.2 


2.4 ±0.2 


X 2 /N 


0.06 


0.21 


0.66 
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Table 7. The same as in Table 6, but for the deuteron. 



Table 7 



parm 


deuteron 


M{ 


Mi 


Mi 




0.022 ±0.001 


0.0060 ±0.0002 


0.0024 ±0.0001 


u n 


0.047 ±0.004 


0.032 ±0.003 


0.020 ±0.002 


/ n 


4.4 ±0.2 


4.2 ±0.2 


4.1 ±0.3 




-0.043 ± 0.003 


-0.030 ±0.003 


-0.019 ±0.002 


In 


2.9 ±0.2 


2.6 ±0.2 


2.4 ±0.3 


X 2 /N 


0.10 


0.37 


1.0 
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FIGURES 




Figure 1. Reduced response functions of the deuteron versus the squared invariant mass W 2 = 
M 2 + Q 2 - (1/x - 1), measured at SLAC @ for Q 2 = 1.75 {GeV/c) 2 (a) and Q 2 = 3.25 (GeV/c) 2 
(b). The open dots and triangles correspond to the (reduced) longitudinal Rl and transverse Rt 
responses, respectively. The solid and dashed lines are the results obtained using the procedure of 
Ref. pi. The dipole form factor G D {Q 2 ) is explicitly given by G D = 1/(1 + Q 2 /0.71) 2 . 
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Figure 2. Pseudo-data for the (Natchmann) transverse moments MT(Q 2 ) (Eq. (|l3|)) of the 
proton versus Q 2 for n = 2, 4, 6 and 8. Open and full dots correspond to the inelastic and elastic 
contributions, respectively, while the open squares are the total moments. The elastic contribution 



is evaluated via Eq. (21) assuming a dipole form for the form factors. 
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Figure 3. The same as in Fig. 2, but for the deuteron. 
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Figure 4. Pseudo-data for the (Natchmann) longitudinal moments M^(Q 2 ) (Eq. fll4|)) of the 
proton versus Q 2 for n = 2, 4, 6 and 8. Open and full dots correspond to the inelastic and elastic 
contributions, respectively, while the open squares are the total moments. The elastic contribution 



is evaluated via Eq. (22) assuming a dipole form for the form factors. 
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Figure 5. The same as in Fig. 4, but for the deuteron. 
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PROTON 





Figure 6. Twist analysis of the proton transverse moments M^(Q 2 ) (Eq. (|l3|)) for n = 2,4,6 and 
8. The open squares represent the pseudo-data of Fig. 2. The thick solid lines are the result of 
Eqs. (|23l) and ( p3|) fitted by the least-x 2 procedure to the pseudo-data of Fig. 2. The thick dashed, 
dotted and dot-dashed lines correspond to the contributions of the twist-2, twist-4 and twist-6, 
respectively. The thin long-dashed lines are the predictions of the twist-2 term obtained using the 
parton distributions of Ref. 32]. 
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Figure 7. The same as in Fig. 6, but for the deuteron. In (a) only the twist-2 term is reported 
(see text). In (b), (c) and (d) the error bars on the different twist contributions correspond to the 
total uncertainties generated by the errors of the parameters appearing in Table 3. 
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Figure 8. The leading-twist moments /m^(Q 2 ) (see Eqs. versus Q 2 for n = 2,4, 6. Open dots: 

twist-2 extracted from our analysis of the proton and deuteron transverse pseudo-data (see also 
Tables 1-3); the errors represent the uncertainty of the fitting procedure corresponding to one-unit 
increment of the x 1 /N variable. Dot-dashed lines: GRV prediction at NLO |p^| . Dashed lines: 
the same as the dot-dashed lines, but including the empirical nuclear correction of Ref. |^3|| . Solid 
lines: the same as the dashed lines, but including the enhancement of the ci-quark distribution at 
large x of Ref. [lilj , given explicitly by d(x) = d(x) + 0.1x(l + x)u(x). 

33 



PROTON 



O 




Cf (GeV/c) 



10 

2 



w q 0.004 




Cf (GeV/c) 



1 o 

2 



w q 0.002 




0.001 



0.000 



■0.001 u 



1 



Cf (GeV/c) 



10 

2 



o 



0.0015 



0.0010 



0.0005 



0.0000 



■0.0005 




Cf (GeV/c) 



10 

2 



Figure 9. Twist analysis of the proton longitudinal moments M^{Q 2 ) (Eq. (|i~4|)) for n = 2,4, 6 
and 8, adopting the Ji2-renormalon model of Ref. fl4]| . The open squares represent the pseudo- 
data of Fig. 7. The solid lines are the results of Eq. (Q) fitted by the least-^ 2 procedure to the 
pseudo-data. The dashed lines are the twist-2 contribution given by Eq. @ adopting the GRV 
set of parton distributions |3^] at NLO. The dotted and dot-dashed lines correspond to the twist-4 
and twist-6 terms, respectively. 
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Figure 10. The same as in Fig. 9, but for the deuteron. 
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Figure 11. Twist analysis of the proton transverse moments M^(Q 2 ) (Eq. (|l3|)) for n = 2,4,6 
and 8, including both phenomenological higher twists originating from multiparton correlations 



and the power corrections arising from the Ii?-renormalon model of Ref. [14| (see text). The 
values of the Ii?-renormalon parameters are fixed at A^ R = —0.132 GeV 2 and A\ R = 0.009 GeV 4 , 
as resulting from the analysis of the longitudinal moments. The solid lines are the results of Eq. 



(34) fitted by the least-x 2 procedure to the pseudo-data (open squares). The dashed lines are the 
twist-2 contribution, while the dotted and dot-dashed lines correspond to the 7\R-renormalon and 
multiparton correlation contributions, respectively. 
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Figure 12. The same as in Fig. 11, but for the deuteron. In (a) only the twist- 2 term and the 
/i?-renormalon contribution are reported (see text). 
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Figure 13. The mass scale (see Eq. (|37|)) of the twist-4 and twist-6 terms resulting from our 
final analysis of the transverse pseudo-data (see Tables 6 and 7). Open and full markers correspond 
to the proton and deuteron case, while dots and squares are our results for the twist-4 and twist-6, 
respectively (see text). 
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Figure 14. The sum of the NLO twist-2 (Eq. (0)) and the LR-renormalon (Eq. (|27|)) contributions 
to the transverse moments M^(Q 2 ) with n > 4, as extracted from our final analysis of the transverse 
pseudo-data (see text), versus Q 2 . Thin and thick lines are the results obtained at a s (M|) = 0.113 
and 0.118, respectively, while the solid, dashed and dot-dashed lines correspond to the transverse 
moments of order n = 4, 6 and 8, respectively. 
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